Relaxation and reconstruction on (111) surfaces of Au, Pt, and Cu 
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We have theoretically studied the stability and reconstruction of (111) surfaces of Au, Pt, and Cu. 
We have calculated the surface energy, surface stress, interatomic force constants, and other relevant 
quantities by ab initio electronic structure calculations using the density functional theory (DFT), 
in a slab geometry with periodic boundary conditions. We have estimated the stability towards 
a quasi-one-dimensional reconstruction by using the calculated quantities as parameters in a one- 
dimensional Frenkel-Kontorova model. On all surfaces we have found an intrinsic tensile stress. This 
stress is large enough on Au and Pt surfaces to lead to a reconstruction in which a denser surface 
layer is formed, in agreement with experiment. The experimentally observed differences between the 
dense reconstruction pattern on Au(lll) and a sparse structure of stripes on Pt(lll) are attributed 
to the details of the interaction potential between the first layer of atoms and the substrate. 
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I. INTRODUCTION 

The properties of close-packed noble and transition 
metal surfaces have been extensively studied in recent 
years. The research has made great advances since the in- 
troduction of new experimental techniques, in particular 
scanning tunneling microscopy (STM) , and the improve- 
ment of first-principle computational methods based on 
the density functional approach. These surfaces show a 
wide variety of behavior with respect to reconstruction, 
preferred site and strength of chemisorption of reactive 
species, adsorbate diffusion, etc. 

Owing to the abrupt change of the electronic struc- 
ture on metal surfaces the first layer of atoms may have 
rather different properties from the bulk. Often, a large 
intrinsic tensile surface stress appears, which is the driv- 
ing force leading to compressive reconstruction on many 
close-packed metal surfaces. The stability of a particular 
surface is the result of the interplay of several physical 
quantities characteristic of surface, such as surface en- 
ergy, surface stress, interatomic force constants, etc. The 
reconstruction is more likely to occur the larger the in- 
trinsic stress and the smaller the energy for sliding the 
atoms of the first layers into positions out of registry with 
the substrate. 

In this paper we consider (111) surfaces of copper, plat- 
inum, and gold. The Au(lll) surface reconstructs form- 
ing a uniaxially compressed layer, which can be seen in 
STM experiments^^ The Pt(lll) reconstructs only at 
high temperatures or in the presence of saturated Pt 
vapourj&i and Cu(lll) does not reconstruct. We report 
on first-principle numerical calculations of the proper- 
ties of these surfaces using the density functional theory 
(DFT). In Sec.|n]we first describe the generalities of the 
computation procedure. We then calculate the important 
physical quantities, such as the effective force constants 
in the first layer, the surface energies, and the intrinsic 
surface stress. In Sec. lIIII we use these quantities in order 



to estimate the stability of the first layer of atoms with 
respect to a quasi-one-dimensional reconstruction, which 
can be treated theoretically by the Frenkel-Kontorova 
model. In the Sec. |W]we discuss the results. 



II. CALCULATIONS OF SURFACE 
PROPERTIES 

The reconstructed phase of some close-packed metal 
surfaces may have a periodicity which involves many sur- 
face atoms. For example, the reconstruction of Au(lll), 
one of the surfaces considered in this paper, involves a 
uniaxially compressed first atomic layer with a period of 
around 70 A. Full first-principle calculations with such 
large supercells are not feasible. One instead evaluates 
the stability by calculating the relevant quantities of the 
unreconstructed surface, and uses them in phenomeno- 
logical models of reconstruction. 

There have been a number of papers in which the 
stability of (111) surfaces of noble metals was evalu- 
ated using the Frenkel-Kontorova (FK) model^^*^ of 
uniaxially compressive reconstruction. Ravelo and El- 
Batanouny (Ref.0, see also references therein) have con- 
structed effective potentials between surface atoms and 
performed molecular-dynamics simulations of the recon- 
structed phase. Mansfiled and Needs^ have derived con- 
ditions of stability of (100) and (111) surfaces towards a 
compressive reconstruction in the FK model, and eval- 
uated it for several metals. Apparently, their values of 
the parameters are not good enough, as they have found 
no reconstruction of Pt(lll) and Au(lll). In particu- 
lar, the strength of the interatomic potentials seems to 
be overestimated. Takeuchi et al*2i have constructed a 
two-dimensional potential between first- and second-layer 
atoms of a Au(lll) surface using the results of first- 
principle calculations. They have calculated the recon- 
struction pattern using simulation techniques in a two- 
dimensional FK model. Using the two-dimensional po- 
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tential in a molecular-dynamics simulation, Narasimhan 
and Vanderbilt — have concluded that the herringbone 
pattern of the reconstruction stripes observed on Au(lll) 
is favored by the long-range elastic interactions mediated 
by the substrate. Recently, Narasimhan and Pushpa 25 
have studied a two-dimensional FK model of the Pt(lll) 
surface for which the parameters have been obtained from 
ab initio calculations. They have obtained simulated 
STM images of the reconstruction stripes recombining 
and intersecting in various ways, similarly to the struc- 
tures observed in STM experiments. 

The aim of this paper is to re-evaluate with highest 
accuracy the properties of the surfaces which determine 
the stability towards reconstruction. In this section we 
report on such calculations for the (111) surfaces of Au, 
Pt, and Cu in the density functional theory (DFT) ap- 
proach. In the following subsections we describe, in or- 
der, the generalities of the numerical calculations, the cal- 
culation of surface energies, surface stress, surface spring 
constants, and the calculation of energy required to slide 
the complete first layer to various positions away from the 
most stable one. In the next section we use these results 
as input parameters in the (effectively one-dimensional) 
Frenkel-Kontorova model of uniaxially compressive re- 
construction. 



A. DFT calculations 

We have performed first-principle calculations of the 
electronic structure of (111) surfaces within the density 
functional theory approach, using the dacapo program. 12 
We have used the provided ultrasoft pseudopotentials for 
the Perdew-Wang exchange-correlation functional PW91 
and the generalized gradient approximation (GGA). 

We have first made calculations for bulk metal, in or- 
der to determine the lattice constant at which the energy 
per atom is minimum, using the primitive fee unit cell. 
The lattice constant found often differs slightly from the 
experimental value: for example, our value is 4.00 A com- 
pared with the experimental value 3.92 A for Pt, 4.18 A 
compared with the experimental value 4.08 A for Au, and 
3.66 A compared with 3.61 A for Cu. We have followed 
the common practice using the computed value of the 
lattice constant in susequent calculations, which ensures 
that no spurious bulk stresses appear. 

The surfaces were described by a slab of five or more 
hexagonal layers of atoms. Since the calculation assumes 
periodicity in all three dimensions, the metallic slabs were 
separated by typically five layers of vacuum. In addition 
to ideal surfaces, we calculated surfaces perturbed in var- 
ious ways in order to deduce quantities like surface stress 
and surface energy, which can be used to estimate the 
stability of the surface. More details are given in the 
respective subsections. 

In the case of a clean surface, the unit cell in the di- 
rections along the surface plane (x — y) consisted of one 
atom, and 56 /c-points in the two-dimensional first Bril- 



louin zone were used. In some calculations described in 
the following subsections, e.g., when alternate rows of 
first-layer atoms were displaced in opposite directions, 
in order to probe the restoring forces within the first 
layer, unit cells with four atoms in each plane were 
used, and the number of fc-points was reduced to 18. In 
the ^-direction, perpendicular to the layers, only k = 
was considered, consistent with the assumption that the 
slabs, which repeat periodically because of the compu- 
tational method, did not couple to each other. In most 
calculations, the two bottom layers were kept fixed at 
the bulk separation. We used an energy cutoff for the 
plane-wave basis set of 340 eV, and the electronic oc- 
cupation was smeared by a pseudothermal distribution 
of T = 0.2 eV. We performed some checks with a lower 
value, T = 0.1 eV, and found that the changes of the 
calculated quantities were irrelevant. 

Next, we calculated the structure of clean (111) sur- 
faces. The relaxation of the surface layers from truncated 
bulk positions was found to be rather small for the three 
metals considered. For gold, the first and second layers 
relax less than 0.3% of the interlayer distance in the (111) 
direction, with an energy gain of 1 meV, which is irrel- 
evant considering the accuracy of the calculation. The 
first platinum layer relaxes outwards by about 1%, with 
an energy gain of 2.5 meV. This unusual expansion of the 
interlayer distance has been observed experimentally. 13 
On copper, the first three layers move inwards by less 
than 1%, and the energy gain is 4 meV. 



B. Surface energy 

We obtained the surface energy per atom as the dif- 
ference of the energy of the bulk and of a slab. In order 
to minimize the influence of the computational details, 
we calculated the energy of the bulk using a six layer 
slab (i.e., twice the abc stacking of the (111) fee direc- 
tion) with no vacuum layers. After that, we changed the 
configuration by introducing four layers of vacuum which 
created two (111) surfaces, and calculated the energy, al- 
lowing both surfaces to relax. The surface energies 7, 
one half of the difference of the energies obtained in the 
two calculations, are given in the GGA column in Ta- 
ble HJ These values are consistently smaller by about 
one third than those calculated recently by Vitos et al^ 
and Galanakis et a\M^ and various experimental val- 
ues reported therein. One possible source of the discrep- 
ancy are the different pseudopotentials and exchange- 
correlation functionals used. We therefore also report the 
surface energies calculated using the LDA functional for 
the electronic densities obtained in the GGA calculations, 
in the column labeled as LDA in TablcQ] These electronic 
densities are, of course, non-self-consistent with respect 
to the LDA functional. They are always larger than the 
GGA results, in better agreement with Refs. ll4TIT5llT^ . 
Another possible source of discrepancy is the small thick- 
ness of the slab in our calculation. In the bulk calculation 
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TABLE I: Surface energy 7, surface stress r, and force con- 
stants between nearest neighbors fco, for (111) surfaces of Au, 
Pt and Cu. The two values of the surface energy 7 have been 
obtained from a fully self-consistent DFT GGA calculation, 
and from applying the LDA functional to the same electronic 
density, see text. The tree values of the force constant fco have 
been calculated from the bulk modulus B, (feos), and from 
the forces obtained in DFT calculations in which the first- 
layer atoms were slightly displaced "compressively" (hoc) or 
"laterally" (k os ), see Fig. [I] 





7 (eV/A 2 ) 


t (eV/A 2 ) 




k (ev/A 2 ; 


1 




GGA 


LDA 




koB 


koc 


kos 


Au 


0.050 


0.080 


0.15 


1.83 


1.41 


1.18 


Pt 


0.084 


0.123 


0.34 


3.12 


2.23 


2.55 


Cu 


0.080 


0.112 


0.11 


2.04 


1.73 


2.36 



bulk term, as discussed earlier. The diagonal element of 
the surface stress tensor r, which is the relevant quantity 
for the surface reconstruction, is 



where AE is the change of energy when the lattice con- 
stant is varied (after subtracting the quadratic bulk term) 
and AA is the associated change of the surface area. The 
numerical values are given in Table [n] We have found 
that in this case there is almost no change if the non-self- 
consistent LDA values of the energy are used instead of 
the GGA values. 



D. Surface spring constants 



(i.e., 6 layers without vacuum) we used only one fc-point 
in the z-direction, in order to keep the similarity with the 
surface calculations, and consequently the long-distance 
contributions in the ^-direction might not be correctly 
taken into account. This can be investigated by using 
elementary cells with more layers, i.e., repeating the abc 
structure three or more times, but we have not done such 
checks. Also, we used the equilibrium lattice constants 
obtained self-consistently in our GGA calculations which 
were larger than the experimental values and those ob- 
tained by other authors from LDA calculations. 

In the following section, when discussing the stability 
of the surfaces, we shall use the LDA results, which seem 
to agree well with the best values reported in the litera- 
ture. 



C. Surface stress 

The surface stress can be found by considering the 
change of energy when the lattice constant in the x and y 
direction (i.e., in the surface plane) is varied. In the bulk 
calculations the leading correction to the energy when the 
lattice constant is varied around the equillibrium value 
is quadratic by construction, since we used the lattice 
constant corresponding to the energy minimum. Owing 
to the reduced coordination, the optimum interatomic 
distance in the surface layer is smaller than the bulk lat- 
tice constant, as already mentioned, causing an intrinsic 
tensile surface stress ^Li£ Therefore, the surface contibu- 
tion to the energy has a leading linear term if the lattice 
constant in the x and y direction is varied. 

In order to find the intrinsic surface stress, we cal- 
culated the energy of slabs consisting of six layers, as in 
the preceding subsection, both unperturbed and with the 
unit cell compressed or expanded along the coordinates 
x and y (i.e., in the directions in the surface plane) for 
around 0.1%, not allowing any additional relaxation in 
the z direction. We also calculated similar configurations 
without vacuum layers, in order to subtract the quadratic 



The simplest atomic-scale model of the lattice dynam- 
ics of fee metals is to assume central harmonic forces 
between nearest neighbor atoms. In this model there 
is a universal scaling of phonon spectra and different ele- 
ments of the elasticity tensor, which is obviously an over- 
simplification. Nevertheless, this model is sufficient for 
the purpose of this work. The force constant kos of the 
nearest neighbor bond in the bulk can be found from the 
bulk modulus^ B = 1/3(C U + 2C 12 ): 



where a is the lattice constant of the conventional fee 
unit cell. 

However, the effective force constant between the 
atoms in the first surface layer can differ from the value in 
the bulk, since the lower coordination can substantially 
alter electronic properties, and an ab-initio DFT calcu- 
lation is necessary. (Note that this force constant also 
contributes to the quadratic energy term in the preced- 
ing subsection, where we did not consider the possibility 
that it was modified within the first layer. However, the 
possible error which this introduces in the calculation of 
the surface stress is negligible.) 

We made slab calculations (with relaxation turned off) 
in which rows of atoms in the first layer were displaced 
by about 0.05 A either in the direction of the rows or 
perpendicularly to it, as shown in Fig. ^ Unlike other 
calculations in this paper, here we were not interested in 
the total energy, but rather in the restoring forces ap- 
pearing owing to the nonequlibrium configurations. We 
analyzed the results assuming central harmonic interac- 
tions ("springs") between nearest neighbors, neglecting 
the coupling to the second layer. Taking into account the 
number of nearest neighbor bonds and the angles, we cal- 
culated the force constants kos- The results for both the 
bulk values obtained from Eq. @ and the surface DFT 
calculations are shown in Table |U We note once again 
that this approach does not attempt to give a complete 
description of the lattice dynamical properties of the first 
surface layer. In general, even if the interactions with 
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TABLE II: The parameters of the Frenkel-Kontorova model of uniaxial compressive reconstruction. A is the area per surface 
atom, a is the distance between the fee and bec hollow site on the (111) surface. The other parameters were determined from 
the DFT calculations as described in the text. These are the surface stress r, the non-self-consistent LDA value of the surface 
energy jlda, the force constant \i — 3/2fcoc an d the- average amplitude of the corrugated potential of the second layer W. 
The quantities a and (3 denned in Eq. £IJ determine the stability of the surface in the Frenkel-Kontorova model. 





A (A 2 ) 


a (A) 


r (eV/A 2 ) 


Ilda (eV/A 2 ) 


H (eV/A 2 ) 


W (eV) 


a 


P 


Reconstruction 


Au 


7.59 


3.41 


0.15 


0.080 


2.115 


0.038 


22.90 


-23.97 


Compressive 


Ft 


6.94 


3.27 


0.34 


0.123 


3.345 


0.061 


21.80 


-37.58 


Compressive 


Cu 


5.79 


2.99 


0.11 


0.112 


2.595 


0.056 


18.32 


-3.58 


No 




fa) (b) 

FIG. 1: The calculation of the force constants between the 
atoms in the surface layer of (111) surface. The atoms in 
the top layer of the slab were displaced by around 0.1% of 
the interatomic distance, and the restoring forces were calcu- 
lated. Assuming that only central nearest-neighbor restoring 
forces exist, the force constants are hoc = (1/6)(F/Sy) for 
"compressive" displacements (a), and fcos = (1/2)(F/Sx) for 
"lateral" displacements (b). 



subsurface atoms are not taken into account, a set of force 
constants between second and further neighbours and/or 
angular force constants would be necessary in order to re- 
produce accurately the lattice dynamics, i.e., the phonon 
spectra. The present approach is a simplification to be 
used only in an estimation of the stability of the surface 
layer. One may further object that the reconstruction 
has a large wavelength, while the displacements shown 
in Fig. ^ correspond mainly to short-wavelength modes 
around the edge of the Brillouin zone, and the estimated 
ks need not be the same. In our opinion, the differences 
are small. 



E. Potential between the top layer and the bulk 

When reconstruction occurs, some first-layer atoms are 
displaced to energetically less favorable positions with re- 
spect to the second layer. The last quantity which we 
need for an estimate of the stability is the amount of en- 
ergy lost by the atoms when displaced to various nonopti- 
mal positions. Since the periodicity of the reconstruction 
pattern is large compared with the substrate periodicity 
(e.g., by a factor of around 22 on Au(lll)), the posi- 
tion of each subsequent atom along the reconstruction 
direction with respect to the second layer changes only 




FIG. 2: Various positions of the first-layer atoms (small cir- 
cles) on a (111) surface, relative to the second layer (large 
circles) and the third layer (grey). In the unreconstructed 
phase, all first-layer atoms are in the fee positions (/). Upon 
reconstruction, the atoms are found in various positions along 
the path f — b— h — b — f. 



slightly. A good estimate of the energy involved can be 
obtained by considering structures in which all atoms in 
the first layer are simultaneously displaced by the same 
vector. We performed DFT calculations of such configu- 
rations, keeping the displacement in the x — y plane fixed 
and allowing the atoms to relax in the z direction. 

We first calculated the energy of the regular fee con- 
figuration, denoted by / in figures and tables, which is 
energetically the most favorable. Next, we considered 
the configuration with the first layer atoms in the hep 
hollows of the second layer (h), which is also a local en- 
ergy minimum. We furthermore calculated the energy 
of some other configurations which are not local energy 
minima, keeping the x and y coordinates fixed, so that 
the algorithm for the atomic relaxation ignored the lat- 
eral forces acting on the first-layer atoms. These are (see 
Fig. |2J) the "bridge" position (b) in which the first-layer 
atoms are halfway between fee and hep hollows, the on- 
top configuration (t), the configurations (ft) and (ht) 
halfway between (/) and (t), and (h) and (t), respec- 
tively, etc. (Not all calculations were performed for all 
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FIG. 3: The energy needed in order to move the complete 
first atomic layer away from the optimum fee into other con- 
figurations. The various positions are depicted in Fig. [5] and 
the energies are expressed in eV/atom. 



three metals.) For the quasi-one-dimensional compres- 
sive reconstruction, the relevant configurations are along 
the path / — fb — b — bh — h and back to b and /. (The 
on-top configuration, which is not occupied in the re- 
construction, was calculated with the aim to obtain a 
better insight into the difference between various met- 
als.) In practice, the calculations were performed so that 
initially, all layers except the top one were kept fixed at 
the bulk configuration, the top layer had fixed x and y 
coordinates, and the z coordinate was allowed to relax. 
In a second step, the bottom two layers were kept fixed, 
the intermediate layers were allowed to relax in all di- 
rections, and the top layer was allowed to relax in the 
z direction only. The second step produced only minor 
changes for the symmetrical configurations (/), (h), (t), 
and (b), but the nonsymmetric configurations changed 
significantly, as the atoms in the second layer (and to a 
lesser degree in deeper layers) relaxed laterally under the 
force exerted by the first-layer atoms. 



In Fig. we show the calculated energies, measured 
with respect to the energy of the fee configuration, which 
was the lowest for all three metals. The "bridge" posi- 
tion, midway between fee and hep hollows, is approxi- 
mately a saddle point. All points were calculated for five 
layers (first one fixed in x and y, two free, two fixed), 
and some points were also calculated for six and seven 
layers. It was found that the energies changed by sev- 
eral meV, but the qualitative behavior remained similar, 
and the estimates of stability towards reconstruction in 
the following section are not affected. Generally, the en- 
ergies of symmetric configurations (h, b, t) increase with 
increasing number of layers, while those of nonsymmetric 
configurations, like (fb) or (ht), decrease. This can be at- 
tributed to the fact that in the latter case there are more 
atomic layers which are free to relax laterally, leading to 
a more complete relaxation. 

The first-layer atoms in positions other than the reg- 
ular fee one (/) are also protruding further out. For all 
three metals, the values are similar, namely, around 0.1 A 
for the bridge (b), 0.02 A for the hep (h), and 0.3-0.4 A 
for the on-top position (t). 

The relevant quantity for the stability estimate is the 
average energy W along the path / — b — h — b — f, in 
the approximation of making a Fourier expansion and 
keeping only the two lowest terms: 

W = AE(b) + AE(h)/2. (3) 



III. STABILITY AND RECONSTRUCTION 

In this section we use the calculated quantities to dis- 
cuss the stability of the surfaces of platinum, gold and 
copper. We first introduce the one-dimensional Frenkel- 
Kontorova model, and then discuss, in order, each of the 
considered surfaces, calculating the stability criteria, and 
comparing the results with the known experimental find- 
ings. 



A. Reconstruction of (111) surfaces 

As already said, the electronic structure of closed- 
packed metal surfaces is very different from that in the 
bulk, with an abrupt decrease in density of conduction 
electrons. The consequences may be a change of the 
length of the bonds between surface atoms, occasionally 
accompanied by a reconstruction which involves a change 
of the number of atoms in the first surface layer. 

The reconstruction of the fee (111) surface involves a 
large number of atoms in inequivalent positions with re- 
spect to the underlying layer of atoms. In scanning tun- 
neling microscopy experiments, the reconstruction of the 
Au(lll) surface can be clearly seen as bright stripesji^S 
owing to the increased height of atoms which are out of 
registry with the second layer. On a larger scale, the 
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stripes form a herringbone pattern, as the quasi-one- 
dimensional reconstruction proceeds along three equiv- 
alent directions on the surface. A typical size of the re- 
constructed Au(lll) supercell is 70 A x 280 A. The re- 
construction of the Pt(lll) surface is similar, but occurs 
only at high temperatures^ or in the presence of sat- 
urated Pt vapouri£i£ STM micrographs show that the 
stripes do not form continuous patterns, but are instead 
well separated, with various types of intersections^ The 
stripes in both systems consist of quasi-one-dimensional 
compressive reconstruction^ and can be treated as 
solitonic solutions of the Frenkel-Kontorova modeli2ii22 
The large-scale two-dimensional structure of the recon- 
structed phase depends upon details of the interactions, 
and will not be discussed here. 



B. The Frenkel-Kontorova model 

The one-dimensional Frenkel-Kontorova modelS*^^ 
consists of a linear chain of atoms, subject to two compet- 
ing interactions with different intrinsic periodicities. The 
atoms interact via a nearest-neighbor harmonic coupling 
with a preferred lattice constant b. In addition, there is 
an external periodic potential (i.e., the potential of the 
second atomic layer) of a periodicity a. Depending on the 
strength of the external potential and the magnitude of 
the spring constant [a, various stable solutions are possi- 
ble. At small "pressures" (i.e., small difference of a and b, 
small /i) the external periodic potential dominates and all 
atoms are in potential minima, forming a commensurate 
phase. As the "pressure" increases, the natural periodic- 
ity of the atomic chain becomes more important. If, say, 
b < a, solitons appear in which the atoms are closer to 
each other, thus gaining energy from the harmonic inter- 
action but paying the cost of increased potential energy 
in the external potential. Finally, when the external po- 
tential is weak compared with the elastic energy needed 
to stretch the atomic chain, the atoms follow the period- 
icity b, forming an incommensurate phase. 

Mansfield and Needs^ have applied the Frenkel- 
Kontorova model to the reconstruction of (111) surfaces 
of the fee metals. They have found that the relevant 
quantity is 

= A(j - 4t/3) = Wp 

where the parameters of the Frenkel-Kontorova model 
have already been expressed in terms of the physical 
properties of real (111) surfaces. Thus, 7 is the surface 
energy, r is the surface stress, fj, is 3/2 of the surface force 
constant, and W is the average potential energy. (The 
factors 4/3 in the stress term and 3/2 in the definition 
of \x appear because the path / — h — / is not straight.) 
a is the distance from the fee hollow site to the nearest 
hep hollow site on the (111) surface, and A is the sur- 
face area per atom. These quantities are the same (or 



closely related to) those calculated in the preceding sec- 
tion for real metal surfaces. There is no reconstruction 
for \P\ < 1, the reconstruction is compressive for P < — 1 
and expansive for P > 1. The quasicontinuum approach 
is valid if the magnitudes of the dimensionless quantities 
a and (3 defined in Eq. J3J are large, say larger than 5. 

In the following subsections we apply this analysis to 
the (111) surfaces of gold, platinum, and copper. The 
relevant quantities derived in Section [H] are summarized 
in Table ITU and the resulting a and (3 are given. 

C. Au(lll) 

In our calculations the first layer of an ideal Au(lll) 
surface relaxes outwards by about 0.005 A, less than 
0.3%, which is at the limit of accuracy of the calculations. 
The energy per surface atom of the hep configuration is 
higher than that of the ideal (fee) structure by around 
10 meV per atom, and that of the bridge by 33 meV 
(Fig. |3J) . This is in excellent agreement with the calcula- 
tions of Takeuchi et aliSI and Galanakis et al.fi^*i& who 
used a density functional approach with a mixed basis 
set. These energies are the smallest of all metals consid- 
ered here. The surface energy is also small, and owing to 
a moderately large tensile surface stress we obtain that 
(3 is slightly larger than a, indicating that the surface 
reconstructs. Experimentally, the surface appears par- 
ticularly prone to reconstruction, and a dense pattern of 
stripes of 23 x -\/3 reconstruction is observed at room 
temperatures with STMi 1 ^ 21 ? 

D. Pt(lll) 

Our calculations show that the first atomic layer on 
Pt(lll) relaxes outwards by 0.023 A or 1%. This some- 
what exceptional behavior has also been found in other 
calculations and confirmed experimentally^ It has been 
attributed to an outward pointing electrostatic force on 
the positively charged atoms of the first layer owing to 
the kind of spill-out of electrons away from the surfaced 
The energy of the bridge configuration is around 46 meV 
larger than that of the regular fee surface. The energy 
of the hep configuration is only slightly smaller than the 
bridge (even less so in calculations with 6 and 7 layers), 
which is different from the other surfaces considered. The 
force constant fc c calculated assuming a "longitudinal" 
displacement of rows of first-layer atoms is reduced by 
about 30% compared with that derived from the bulk 
compressibility. The intrinsic tensile surface stress is 
large, more than twice larger than in the other two met- 
als considered. Although the surface energy is also quite 
large, the quantity (3 is much larger than a, predicting a 
strong tendency to reconstruct. 

Experimentally, the reconstruction of Pt(lll) is ob- 
served only at high temperatures or in the presence of 
saturated platinum vapor The nature of the recon- 
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struction is similar to that of Au(lll), but the stripes 
of compressive "solitons" do not form a dense pattern. 
Instead, they are sparse, and occasionally intersect in 
several distinct waysAS^ 

E. Cu(lll) 

Unlike the other two metals considered, we find that 
the first layer on an ideal (fee) Cu(lll) surface relaxes 
inwards by 0.014 A. The contraction of the interlayer dis- 
tance is usual for close packed surfaces of many transi- 
tion metals. The hep configuration is a clear local energy 
minimum, and the energy of the on-top configuration is 
found to be exceptionally large at about 250 meV. The 
intrinsic surface stress is again tensile, but rather small. 
The calculated (3 is clearly smaller than a, indicating that 
the surface does not reconstruct. Indeed, no reconstruc- 
truction has been observed experimentally. 

F. Discussion 

The use of the stability condition of (111) surfaces de- 
rived in the one-dimensional Frenkel-Kontorova model, 
Eq. I|4"|). gives correct qualitative answers, namely, that 
the surfaces of gold and platinum reconstruct compres- 
sively and the surface of copper does not. From a quan- 
titative point of view, the agreement between theory and 
experiment is less satisfactory. We find that the condi- 
tion for reconstruction, \/3\ > a, is barely satisfied for the 
Au(lll) surfaces, while in experiments a dense pattern of 
reconstruction stripes is always found on this surface. On 
the other hand, we obtain that the condition of Eq. (@J 
is amply fulfilled for the Pt(lll) surface, while experi- 
mentally, reconstruction is observed only under favorable 
thermodynamic conditions, at high temperatures, or in 
the presence of Pt vapor. It has already been noticed that 
calculations predict only a marginal tendency for recon- 
struction on Au(lll), at variance with experiment. It has 
been suggested that the ordering of the reconstruction 
stripes into a herringbone pattern is indeed an essential 
contributing factor to the stability of the reconstructed 
phase. 20 ' 22 The theoretical overestimate of the tendency 
to reconstruct in the case of Pt(lll) may be due to the 
large value of the surface energy, which makes it unfa- 
vorable to bring additional atoms to the surface once the 
surface stress has been reduced by the formation of a few 
sparse reconstruction stripes. For a more complete treat- 
ment, a two-dimensional Frenkel-Kontorova model with 
a more realistic potential, calculated in a region around 
the minima, should be used. 

The approach used in this paper can also give relevant 
parameters for other physical properties. We have cal- 
culated relaxations around adatoms of the same species 
on the surfaces considered in this paper, and found that 



the force-relaxation ratio is consistent with the model of 
elastic constants calculated in Section |2] Surprisingly, 
in our calculations (one adatom per four surface atoms) 
the adsorption into a hep hollow site was slightly more 
favorable than into a fee hollow site for all three metals 
considered. We have not found any conclusive experi- 
mental data about this property. 

The approach can also be applied to relaxation around 
defects and chemisorbatcs of different species, energetics 
of self- and heterodiffusion, and other properties of (111) 
surfaces. 



IV. CONCLUSIONS 

Using the density functional theory approach we have 
calculated key properties of (111) surfaces of several met- 
als, and evaluated the stability towards reconstruction 
into a uniaxially compressed reconstructed layer, using 
a one-dimensional quasicontinuum Frenkel-Kontorova 
model. The intrinsic surface stress is tensile for all sur- 
faces, i.e., the atoms in the first layer prefer a denser 
packing than dictated by the potential of the second 
layer, which reflects the bulk periodicity. On the opposite 
side of the energy balance is the energy required to bring 
an extra atom into the surface layer, and the energy lost 
owing to the fact that surface atoms in a reconstructed 
phase are no longer in the minima of the potential of the 
second layer. We have found that Pt(lll) reconstructs 
owing to a large intrinsic tensile surface stress. The stress 
on Au(lll) is considerably smaller, but the surface en- 
ergy is also smaller and the reconstruction criterion is 
satisfied. The intrinsic surface stress in Cu(lll) is some- 
what smaller than that in Au(lll), but other quantities 
are unfavorable, and the surface does not reconstruct, in 
agreement with experimental observations. Experiments 
show that Au(lll) has the largest tendency to recon- 
struct, and dense reconstruction patterns are observed 
using STM, while Pt(lll) reconstructs only under fa- 
vorable conditions, and only a small fraction of the sur- 
face is reconstructed. This is the opposite from what our 
reconstruction criteria suggest, where Pt(lll) is much 
deeper in the reconstruction regime. The reason may 
be further stabilization of the reconstruction of the gold 
surface by the formation of a two-dimensional herring- 
bone structure)22i22i2ii which is not included in the one- 
dimensional model used in the derivation. 
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